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Abstract We calculated the nucleon self-energies in iso-spin asymmetric nuclear matter and obtained the nuclear 


symmetry energy by taking difference of these of neutron and proton. We find that the scalar (vector) self-energy part 


gives a negative (positive) contribution to the nuclear symmetry energy, consistent with the result from relativistic 


mean-field theories. Also, we found exact four-quark operator product expansion for nucleon sum rule. Among them, 


twist-4 matrix elements which can be extracted from deep inelastic scattering experiment constitute an essential part in 


the origin of the nuclear symmetry energy from QCD. Our result also extends early success of QCD sum rule in the 


symmetric nuclear matter to the asymmetric nuclear matter. 
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1 Introduction 


Recently the Korea government decided to build low 
energy rare isotope accelerator, RAON, which will be 
mainly used to study exotic properties of iso-spin 
asymmetric nuclear matter ranging from rare isotope 
to neutron star. To understand the exotic property of 
the asymmetric nuclear matter, one should know what 
cause the difference with the symmetric nuclear 
matter. 

In Dirac phenomenology, strong vector 
repulsion and scalar attraction in optical potential is 
needed to explain the phenomena in nucleon-nucleus 
scattering and these tendency naturally comes from 
successive relativistic mean field type model 
calculations. 

Almost 20 years ago, this tendency had been 
confirmed in QCD degrees of freedom, using finite 
density QCD sum rule, by Cohen, Furnstahl and 
Griegel" “, These studies show that nucleon bulk 
properties can be understood in QCD degree of 
freedom. 

With these two motivations and in a hope to 
express the origin of the nuclear symmetry energy in 
QCD degree of freedom, we have applied QCD sum 


rule to the asymmetric nuclear matter. 


2 Brief description for theoretical method 


2.1 Phenomenological approach 


There was early attempt which describe iso-spin 
asymmetry in a finite nuclei by taking linear potential. 
Following is well known Bethe-Weizsäker formula: 


Eg =a/A-asA3 -ac |Z (Z -1)]A 3 (1) 


—a,I°A+ 6(A,Z), 


where J=(N—Z)/A, the fourth term in Eg accounts for 
the total shifted energy due to neutron number excess. 
When the nuclear symmetry energy is assumed to 
depend on linear density term, the fourth term can be 
easily calculated as 

1 
Ar 


For infinite nuclear matter case, energy 


a,=—|E,(¢,1)-E, (p,1)]. (2) 


density per nucleon of asymmetric (neutron rich) 
nuclear matter can be written as follows: 


E = Eppa +E Pp =E) +E gm (P? +, (3) 
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where E, (E,) is the averaged quasi-nucleon energy 
in asymmetric nuclear matter which can be obtained 
by density integration of single nucleon energy En. If 
one assume linear density dependence of E, (Ep), it is 
easy to get the averaged single nucleon energy. 

By mean-field type approximation which gives 
plausible and clear view for the single nucleon energy 
in the nuclear matter, a single nucleon in the 
asymmetric nuclear matter can be regarded as a 
quasi-nucleon on the asymmetric Fermi-sea with large 
self-energies. From the difference of each nucleon 
self-energy, the nucleon symmetry energy can be 
obtained. 


2.2 QCD sum rule for nucleon in asymmetric 


nuclear matter 


In this study, we followed the formalism in Ref.[1-4] 

and detailed descriptions such as equations and 

formulas are given in Ref.[5]. The brief steps are as 
follows. 

1) QCD sum tule starts with calculating the 
correlation function. The correlation function 
contains all possible quantum states of 
quasi-nucleon in the asymmetric nuclear matter. 

2) The 


calculated by operator product expansion (OPE). 


correlation function can be directly 
The OPE of correlation function can be written as 
multiples of Wilson coefficient and localized 
operator expectation value, the condensate. At 
short distance, the Wilson coefficient can be 
calculated in perturbative diagram calculation. 

3) Also, the correlation function is expected to have 
spectral structure similar as the propagator in 
relativistic mean field theory (RMFT). In this 
phenomenological ansatz, the nucleon 
self-energies are located near the quasi-nucleon 
pole state. 

4) Because we are only interested in the 

self-energies which live near the quasi-nucleon 

pole, the OPEs and phenomenological side 


should be weighted. After the weighting process 


named as Borel transformation, the pole 
contribution in both expression becomes 
enhanced and the continuum contribution 


becomes suppressed. 
5) By equating both weighted expression, one can 


relate the nucleon self-energies with OPE terms, 
in QCD degrees of freedom. 


3 Nuclear symmetry energy from QCD 
sum rules 


To check the physical quantity calculated from QCD 
sum rule, one should set appropriate range of Borel 
mass. It requires following criteria: 

1) The contribution of highest mass dimension 
condensate should not excess 50% of the 
contribution of lowest mass dimension 
condensate. 

2) The continuum contribution of each OPE term 
should not excess 50% of the pole contribution. 

The nuclear symmetry energy which contains 
the condensate up to mass dimension-five is shown in 
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Fig.1 (Color online) Nuclear symmetry energy which contains 
up to mass dimension-five condensates. 

From Fig.1, one can check that the nuclear 
symmetry energy has a value near 40 MeV, consistent 
with previous studies in order of magnitude. 
nuclear 


In following, we present the 


symmetry energy which contains up to mass 
dimension-six condensates. The section is divided into 
the two cases, one contains twist-4 matrix elements 
and the other does not. 

So far, in nucleon OPE, no exact contribution 
of four-quark operators has been known. In this study 
we have found the exact four-quark operators. 
Especially, among them, dimension-six spin-2 
operators (twist-4) can be extracted from DIS 
experiment! and their contribution to the nuclear 


symmetry energy is so important. 
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3.1 Results without twist-4 operator contribution 


If one just account the mass dimension-six spin-1 
condensates, the result is as following figure. 
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Fig.2. (Color online) Nuclear symmetry energy which contain 
up to mass dimension-six condensate (without twist-4). This 
result deals same result presented in Phys Rev C, 87: 015204. 

From Fig.2, one can find the weak vector repulsion 
and scalar attraction consistent with the RMFT result 


which includes iso-vector meson exchange channel”. 


3.2 Results with twist-4 operator contribution 


By including the twist-4 contribution, one can check 
that the mechanism of repulsion-attraction becomes 
more apparent one from Fig.3. This result shows that 
the reaction mechanism in the iso-spin asymmetric 
nuclear matter may be originated from multi-quark 
operators such as twist-4 operators. 
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Fig.3 (Color online) Nuclear symmetry energy which contains 
up to mass dimension-six condensate with twist-4 matrix 
elements. This result deals same result presented in Phys Rev C, 
87: 015204. 


4 Conclusion 


In this study, we calculated the nucleon self-energy in 
asymmetric nuclear matter in QCD degree of freedom 
by using finite density QCD sum rule technique. We 
successfully reproduced the nuclear symmetry energy 
of previous studies in order of magnitude and found 
the exact four quark condensates in the nucleon sum 
rule. From the calculation results, we confirmed 
meaningful fact that the iso-vector meson exchange 
mechanism suggested from many RMFT type studies 
may have origin from multi quark operators such as 
twist-4 operators. The contribution of twist-4 matrix 
element mimics following RMFT results”, 


EJ" = | fa- fa (E) (4) 


The twist-4 matrix elements contain some 
amounts of uncertainties, but the upgrade plan of 
Jefferson Lab. may provide precise twist-4 matrix 
elements and from that, we could get deeper 
understanding of the role of twist-4 operators in the 
nuclear symmetry energy. Still yet the high density 
behavior of the nuclear symmetry energy has not been 
uncovered. Such unknown properties can be studied 
by experimental study, for example by using RAON, 
or by theoretical study. 


References 


1 Cohen T D, Furnstahl R J, Griegel D K. Phys Rev Lett, 
1991, 67: 961-964. 

2 Furnstahl R J, Griegel D K, Cohen T D, Phys Rev C, 1992, 
46: 1507-1527. 

3 Jin X M, Cohen T D, Furnstahl R J, et al. Phys Rev C, 
1993, 47: 2882-2900. 

4 Jin X M, Nielsen M, Cohen T D, et al. Phys Rev C, 1994, 
49: 464-477. 
Jeong K S and Lee S H. Phys Rev C, 2013, 87: 015204. 

6 Choi S, HatsudaT, Koike Y, et al. Phys Lett B, 1993, 312: 
351. 

7 Baran V, Colonna M, Greco V, et al. Phys. Rep, 410: 
335—466. 


050506-3 


